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Abstract 

The cross section of the 16 O(p,p') 16 O(0~, T = 1) scattering was measured at a bombarding en- 
ergy of 295 MeV in the momentum transfer range of l.Ofm -1 < g c . m . < 2.1 fm _1 . The isovector 0~ 
state at E x = 12.8 MeV is clearly separated from its neighboring states owing to the high energy 
resolution of about 30 keV. The cross section data were compared with distorted wave impulse 
approximation (DWIA) calculations employing shell-model wave functions. The observed cross 
sections around q c , m , ~ 1.7 fm^ 1 are significantly larger than predicted by these calculations, sug- 
gesting pionic enhancement as a precursor of pion condensation in nuclei. The data are successfully 
reproduced by DWIA calculations using random phase approximation response functions including 
the A isobar that predict pionic enhancement. 

PACS numbers: 21.60.Jz,25.40.Ep,27.20.+n 
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The search for pionic enhancements in nuclei has a long and interesting history. These 
phenomena can be considered as a precursor of the pion condensation 1| that would be 
realized in neutron stars. Enhancements of the Ml cross section in proton inelastic scattering 
DDDQ and of the ratio Rl/Rt, the spin-longitudinal (pionic) response function Rl to 
the spin-transverse one Rt, in the quasielastic scattering (QES) region 0, Q| were expected 
around a momentum transfer g^m — 1.7 fm -1 . However, the experimental data did not 
reveal any pionic answer the question 

why no pionic enhancements were observed. For example, Bertsch, Frankfurt, and Strikman 
|ll[ suggest the modification of gluon properties in the nucleus that suppresses the pion field. 



Brown et al. 



12j suggest the partial restoration of chiral invariance with density. However, 



we should note that the Ml cross section includes both pionic and non-pionic contributions 
and Rl/Rt is the ratio to the non-pionic Rt- Thus, in these indirect measurements, the 
pionic enhancement might be masked by the contribution from the non-pionic component. 
Recent analyses of the QES data ^| show a pionic enhancement in the spin-longitudinal 
cross section that well represents the Rl, and suggest that the lack of enhancements of 
Rl/Rt is due to the non-pionic component. 

In order to measure the pionic enhancement directly, it is desirable to investigate isovector 
J 71 = 0~, ± — > T excitations because they carry the same quantum numbers as the pion 



and they are free from non-pionic contributions. Orihara et al. 15J] measured the angular 
distribution of the 16 0(p, n) 16 N(0~, 0.12 MeV) reaction at T p = 35 MeV. They reported 
discrepancies between distorted wave Born approximation calculations and their data in the 
range of q c , m , = 1.4-2.0 fm -1 that might be a signature of pionic enhancement. However, 

to the 0~ T = 1 state in 16 at T„ = 65 MeV, such an 



p 



in the proton inelastic scatterin, 
enhancement was not observed 16j. The differences between these (p,n) and (p,p') results 
might indicate contributions from complicated reaction mechanisms at these low incident 
energies. To our knowledge, there are no published experimental data for the 0~, T = 1 
state at intermediate energies of T p > 100 MeV where reaction mechanisms are expected to 
be simple. 

In this Letter, we present the measurement of the cross section for the excitation of the 0~, 
T = 1 state at E x = 12.8 MeV in 16 using inelastic proton scattering at 295 MeV incident 
energy. The results are compared with distorted wave impulse approximation (DWIA) cal- 
culations with shell-model (SM) wave functions. Evidence of a pionic enhancement is clearly 
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observed from a comparison between experimental and theoretical results. The data are also 
compared with DWIA calculations employing random phase approximation (RPA) response 



functions including the A isobar in order to assess the pionic enhancement quantitative 
The measurement was carried out by using the West-South Beam Line (WS-BL) 



III 



and the Grand Raiden (GR) spectrometer at the Research Center for Nuclear Physics, 
Osaka University. The WS-BL provides the beam with lateral and angular dispersions of 
37.1 m and —20.0 rad, respectively, which satisfy the dispersion matching conditions for 
GR. The beam bombarded a windowless and self-supporting ice (H 2 0) target j3] with a 
thickness of 14.1 mg/cm 2 . Protons scattered from the target were momentum analyzed 
by the high-resolution GR spectrometer with a typical resolution of ~30 keV FWHM. The 
beam energy was determined to be 295 ± 1 MeV by using the kinematic energy shift between 
elastic scattering from 1 H and 16 O. The yields of the scattered protons were extracted using 
the peak-shape fitting program allfit 20j • 

The elastic scattering data on 16 O are shown in Fig. El Differential cross sections were 

n 

normalized to the known p + p cross section |21|| by utilizing the data of protons scattered 
from the hydrogen present in the ice target. The data were analyzed using optical model 
potentials (OMPs) generated phenomenologically. The solid curve in Fig.^is the result using 



the global OMP optimized for 16 [22j. The band represents the results by using several 
OMPs parametrized for nuclei from 12 C to 208 Pb with a smooth mass number dependence 
[22) 1 that shows the ambiguity of the OMP for 16 0. The global OMP for 16 reproduces 
the experimental data reasonably well. The systematic uncertainty for the cross section 
normalization is estimated to be less than ~10% from this result. In the following, we will 
use this OMP in DWIA calculations for inelastic scattering. 

Figure El shows the excitation energy spectrum of the 16 0(p, p') scattering at g c . m . = 
1.9 fm . The isovector 0~ state at E x = 12.8 MeV is clearly resolved from the neighboring 
states. The dashed curves represent the fits to the individual peaks while the straight line and 
solid curve represent the background and the sum of the peak fitting, respectively. Narrow 
peaks of 16 were described by a standard hyper-Gaussian line shape, and the peaks with 
intrinsic widths were described as Lorentzian shapes convoluted with a resolution function 
based on the narrow peaks. The positions and widths were taken from Ref. j^J. 

Figure El shows the measured data points and the calculated curves of the cross sections 
of the 0~, T = 1 transition in 16 0(p, p') as a function of the momentum transfer q c . m .- The 
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angular distribution was measured in the range of g c . m . ~ l.Ofm -1 to ~ 2.1 fin -1 starting 
near the second maximum at g c m . ~ 0.9 fm -1 and extending beyond the third maximum at 
q cm ~ 1.7 fm -1 . The error bars of the data points are the fitting uncertainties originating 
from the statistical uncertainties. The shaded areas represent the systematic uncertainties 
including the background subtraction. 

We performed DWIA calculations by using the computer code DWBA98 j^J. The one- 
body density matrix elements (OBDME) for the isovector 0~ transition of 16 0(p, p') were 
obtained from Ref. 



251 ]. This SM calculation was performed in the Os-Op-lsOd-Oflp con- 
figuration space by using phenomenological effective interactions. In the calculation, the 
ground state of 16 was described as a mixture of Ohu; (closed-shell) and 2ftw configurations. 
The single particle radial wave functions were generated by using a Woods-Saxon (WS) 
potential |26j, the depth of which was adjusted to reproduce the separation energies of the 
Opi/2 orbits. The unbound single particle states were assumed to have a very small binding 
energy of 0.01 MeV to simplify the calculations. The AW t-matrix parametrized by Franey 
and Love at 325 MeV was used. The DWIA result is shown as the solid curve in Fig. 01 
The calculation reproduces the data in the lower-g c m . region reasonably well, but they sig- 
nificantly underestimate the data in the higher-g c m . region. Also, the data has a maximum 
at q c , m , ~ 1.7 fm -1 , whereas the maximum of the theoretical curve is slightly higher at q c , m , 
~ 1.8 fm" 1 . 

We investigated the sensitivity of the DWIA calculations to changes of the parameters 
involved. The dash-dotted curve represents the DWIA calculation with a different t-matrix 
parametrized at 270 MeV. The result is systematically larger compared to the calculation 
with the t-matrix at 325 MeV. The dash-dotted curve is, therefore, multiplied by a factor 
of 0.7. The dashed curve denotes the calculation employing a different OBDME with a 
)ure 0p7, 2 lsi/ 2 transition from the Oftw (closed-shell) ground state. Auerbach and Brown 
2fi| suggest that this isovector strength is quenched and spread by a 2f\w admixture. They 
obtained a quenching factor of ~0.64. Thus we have applied this factor as a normalization 
factor to the result. We also performed a DWIA calculation with the radial wave functions 
generated with a harmonic oscillator potential with a size parameter of a = 0.588 fm -1 
|28| . The result is systematically larger compared to the calculation with the WS potential. 
However, their shapes of the angular distribution are very similar to each other. From these 
calculations we found that the shape of the angular distribution is insensitive to changes of 
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the input parameters. Thus it is difficult to understand the discrepancies between experi- 
mental and theoretical results within the framework of the standard DWIA employing SM 
wave functions. Therefore, in the following, we investigate non-locality, RPA correlation, 
and A effects that are not taken into account in these standard calculations. 

The non-locality of the nuclear mean field can be included by introducing a local effective 
mass approximation in the form of 

m*{r) = m N - { WS ^| (m N - m*(0)), (1) 

where is the nucleon mass and /ws( r ) is a WS radial form. The upper panel of Fig. 0] 
shows the m* dependence of the DWIA calculations with the free response function that 

n 

were performed using the computer code CRDW developed by the Ichimura group [29] for 
the analysis of QES data. The 0~ component of the free response is configured as a pure 
0py 2 lsi/2 transition. The DWIA result with m*(0) = is in good agreement with the 
calculation employing the corresponding SM wave function represented by the dashed curve 
in Fig. El Thus we have applied the same normalization factor of 0.64 to the results shown 
in Fig. HJ The angular distribution shifts to lower q c m when decreasing m*(0). As seen in 
the upper panel of Fig. 01 a value of m*{Q) /uin — 0.7 improves the agreement with the data, 
consistent with theoretical estimations [20, till • However there is still a large discrepancy 
between experimental and theoretical results around g c . m . ~ 1.7 fm . 

Next, we discuss the RPA correlation and A effects. We performed DWIA calculations 
with the RPA response functions employing the tt + p + g' model interaction V e s and the 



meson parameters from a Bonn potential which treats the A explicitly 32]. The V e g is the 
sum of the one-7r and one-p exchange interactions, and the Landau-Migdal (LM) interaction 
Vlm specified by the LM parameters, g' NN , q'nai an d 3aai as 



LM 



mi 



A ' V xhyA 9' N A {((ri • T 3 ){tri ■ S 2 ) + h.c.) + (!<-> 2)} (2) 



ml 



+ fr^AA {(Tt ■ T 2 t )(5 1 • Si) + (1 «-> 2)} 



8(ri - r 2 ), 



where cr (t) is the nucleon Pauli spin (isospin) matrix, S (T) is the spin (isospin) transition 
operator that excites N to A, f^NN (/ttVa) is the 7r iViV (7riVA) coupling constant, and m n 
is the pion mass. The middle panel of Fig. 0] shows the q'nn dependence for q'^n = 0.5-0.8 



in 0.1 steps with the fixed g' NA = 0.4 and m*(0)/mN = 0.7. The lower panel shows the g' NA 
dependence for g' NA = 0.2-0.5 in 0.1 steps with the fixed g' NN = 0.7 and m*(0)/mN = 0.7. 



We fixed g' AA = 0.5 |33( since the dependence of the calculated results on this parameter 
is very weak. The calculated angular distributions depend strongly on g' NN in the whole 
g c m . range, whereas a strong g' NA dependence is observed only around g c m . ~ 1.7 fin -1 . The 
most probable choices are g' NN ~ 0.7 and g' NA = 0.3-0.4. As a result, the V e s in the NN 
channel are close to zero whereas that in the iVA channel becomes very attractive . This 
attraction causes the pionic enhancement. 

In conclusion, our high-resolution measurement of 16 0(p, p') 16 O(0~, T = 1) has enabled 
us to search for a pionic enhancement at an intermediate energy of T p = 295 MeV where 
the theoretical DWIA calculations are reliable owing to the simple reaction mechanism. A 
significant enhancement has been observed around q c , m , ~ 1.7 fm -1 compared to standard 
DWIA calculations with SM wave functions. The DWIA analyses employing RPA response 
functions with the A isobar allowed the determination of m*(0)/mN ~ 0.7 and a set of 
the Landau-Migdal parameters of g' NN ~ 0.7 and g' NA = 0.3-0.4, in good agreement with 
estimates based on other experimental data The present direct measurement strongly 
indicates the existence of the pionic enhancement in nuclei that may be attributed to a 
precursor phenomena of the pion condensation. 

We thank the RCNP cyclotron crew for providing a good quality beam. This work was 
supported in part by the Grants-in-Aid for Scientific Research Nos. 12740151, and 14702005 
of the Ministry of Education, Culture, Sports, Science, and Technology of Japan. 
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FIG. 1: The measurement of the cross section for 16 0(p, p) at T p = 295 MeV. The solid curve is 
the theoretical prediction using the global OMP for 16 O. The band represents the ambiguity of 
the prediction as explained in the text. 
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Excitation energy of 16 

FIG. 2: The excitation energy spectrum for 16 0(p,p') at T p = 295 MeV and q c . m . = 1.9 fm _1 . The 
curves show the reproduction of this spectrum with hyper-Gaussian and Lorentzian peaks and a 
continuum. 
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FIG. 3: The measurement of the cross section of 16 0(p,j/) 16 0((T, T = 1) at T p = 295 MeV. The 
shaded areas represent the systematic uncertainties of the data. The solid (dash-dotted) curve 
is the DWIA result with the t-matrix parametrized at 325 (270) MeV employing the SM wave 
function in the Os-Op-lsOd-lpOf model space. The dashed curve denotes the DWIA result with the 
t-matrix at 325 MeV employing the pure 0pi/ 2 ]s±/ 2 SM wave function. 
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FIG. 4: The top panel shows the m*(0)/mN dependence of the DWIA calculations with the free 

response function. The middle panel represents the g' NN dependence of the calculations employing 

the RPA response function with fixed g' NA = 0.4 and m*{Q)/mN = 0.7. The bottom panel denotes 

the g' NA dependence of the calculations employing the RPA response function with fixed g' NN = 

0.7 and m*(0)/mN = 0.7. The curves in each panel represent calculations where the parameters 

in the ranges shown are changed in steps of 0.1. 
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